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Abstract. 

Stellar magnetic fields are produced by a magnetohydrodynamic dynamo mechanism working 
in their interior - which relies on the interaction between plasma flows and magnetic fields. The 
Sun, being a well-observed star, offers an unique opportunity to test theoretical ideas and mod- 
^ ■ els of stellar magnetic field generation. Solar magnetic fields produce sunspots, whose number 

. increases and decreases with a 11 year periodicity - giving rise to what is known as the solar 

, cycle. Dynamo models of the solar cycle seek to understand its origin, variation and evolution 
with time. In this review, I summarize observations of the solar cycle and describe theoretical 
' ideas and dynamo modeling efforts to address its origin. I end with a discussion on the future 
of solar cycle modeling - emphasizing the importance of a close synergy between observational 
data assimilation, kinematic dynamo models and full magnetohydrodynamic models of the solar 
interior. 
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2 1. Introduction to the Solar Cycle 

Sunspots have been observed for many years and systematically so from the early 17th 
^ ^ ^ Century since the invention of the telescope. The number of sunspots on the surface of 
the Sun varies cyclically going through successive maxima and minima vifith an average 
period of 11 years (Fig. 1-top). Given that sunspots are strongly magnetized regions (Hale 
' 1908) with typical field strengths on the order of 1000 Gauss (G), this 11 year solar cycle 
■ is essentially a magnetic cycle. At the beginning of the cycle, sunspots are observed to 
I appear at high latitudes and with the progress of the cycle, more and more sunspots 
ly-j ■ appear at lower and lower latitudes with the cycle ending close to the equator. This 
equatorward migration of the sunspot formation belt takes place in both the hemispheres 
generating the well known solar butterfly diagram of sunspots (Fig. 1-bottom). At the 
same time, the weaker, radial field outside of sunspots are seen to move poleward with 
the progress of the cycle - reaching the poles and starting to reverse the old cycle polar 
^ field around the time of sunspot maxima. The polar field reaches its strongest intensity 
• ^ at sunspot minimum, thereby having a 90° phase lag with the sunspot cycle. The polar 
field reversal occurs with an average period of 11 years as well and is therefore coupled 
^ to the sunspot cycle. 

" " The solar cycle is not strictly periodic. There are variations from one cycle to another 

in both the period and more notably in the amplitude of the cycle. Even the length 
and nature of minima varies from one cycle to another, which was highlighted recently 
in the recent, unusually long solar cycle 23 minimum (see Supplementary Information- 
Fig. 1 from Nandy, Mufioz-Jaramillo & Martens 2011). In the recorded history of sunspot 
observations there is a period between 1645-1715 AD when hardly any sunspots were 
seen on the solar surface (Fig. 1-top). This period is termed as the Maunder minimum 
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Figure 1. Top: The variation in the number of sunspots observed on the Sun with time over the 
last fours centuries. The variation is for the most part, cyclic, ranging typically from 9-14 year 
cycles and with an overall average period of 11 years. The amplitude of the cycle is observed 
to fluctuate significantly from one cycle to another. A notable period is the Maunder minimum 
between 1645-1715 AD, when the sunspot cycle pretty much stopped. Solar activity reconstruc- 
tion based on cosmogenic isotopes indicate that there have been many such episodes of reduced 
or no activity. Bottom: A butterfly diagram (latitude-time plot) of solar surface magnetic field 
showing the evolution of latitudes where sunspots emerge (black regions) and the surface radial 
field (outside of sunspots; scale on right-hand y-axis) with the progress of time. About two and 
half sunspot cycles are depicted here. The butterfly diagram shows that the sunspot emergence 
belt moves equatorward, whereas, the weak surface radial field moves poleward with the progress 
of the solar cycle. The polar field reversal starts at about the time of sunspot maximum, and 
the polar field is the strongest at sunspot minimum. 

and it is thought that the sunspot cycle stopped during this phase but restarted again. 
Solar activity reconstructions based on cosmogenic isotopes indicate that in the past, 
several such episodes have occurred in the Sun (Usoskin, Solanki & Kovaltsov 2007). 
Such grand minima are therefore an integral part of solar activity over long timescales. 
Some evidence based on tree-ring data also suggests that there was a weak cycle - at 
least in the solar open flux and polar fields - during the Maunder minima (Miyahara et 
al. 2004). 
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Figure 2. A magnetogram showing tlic vortical component of the solar magnetic field on the 
Sun's surface. The typical strength of the magnetic field within sunspots is on the order of 1000 
G. White regions denote sunspots which have positive vertical magnetic field and dark regions 
show sunspots with negative vertical magnetic field. Sunspots occur in bipolar pairs which have 
opposite polarities. Their polarity orientation is reversed across the equator, which is evident in 
the above figure. The line joining the two polarities of a sunspot pair are systematically tilted 
relative to the local parallel of latitude and this tilt increases with increasing latitude away from 
the equator. 

The sunspots themselves are observed to emerge in pairs of opposite polarities (see 
Fig. 2). The bipolar sunspot pairs have a systematic tilt with pairs at higher latitudes 
being more tilted than those at low latitudes - following what is known as the Joy's law 
distribution of tilt angles (Hale et al. 1919). It is also seen that the orientation of these 
bipolar sunspot pairs are always reversed across the equator, i.e., if the leading polarity 
spots on the Northern-hemisphere are positive, then the leading polarity spots in the 
Southern- hemisphere are negative. This relative orientation reverses from one 11 year 
cycle to another and therefore the full magnetic cycle constitutes 22 years. 

It is to be noted that similar magnetic cycles - exhibiting a wide range of activity 
behavior including Maunder minima like states are also observed in other solar-like 
stars (Baliunas & Vaughan 1985; Giampapa et al. 2006; Poppenhager et al. 2009). This 
indicates that magnetic cycles are a common and ubiquitous phenomena in stellar physics 
and plausibly have a common underlying mechanism at its origin (see e.g., Saar Bran- 
denburg 1999, Nandy 2004). 

Plasma flows in the Sun are believed to play a crucial role in the solar cycle and 
are now observed routinely. Tracking sunspots was an early indicator of the differential 
rotation of the Sun and it is known that the equator of the Sun rotates faster than 
its poles. Now, with the observational technique of helioseismology we can also observe 
the differential rotation in the solar convection zone (SCZ) and know that the mainly 
latitudinal differential rotation in the SCZ changes over to a radial differential rotation 
at its base and is concentrated in a thin layer known as the tachocline (Charbonneau 
et al. 1999). The differential rotation of the Sun constitutes the ^-component of the 
plasma flow in the solar interior. There is also another component of the plasma flow in 
the Sun's interior in the meridional plane {r-9 plane); this meridional flow of plasma is 
observed to be poleward in the surface and near-surface layers (Gonzalez Hernandez et 
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al. 2006) and is believed to pervade the SCZ (Giles et al. 1997), with an equatorward 
counterflow theoretically expected at or below the base of the SCZ (Nandy & Choudhuri 
2002). Convective turbulence is also expected to be a feature of the plasma in the Sun's 
interior and visible manifestations of super-granular cells in the surface are an indicator 
of the underlying convective turbulent motions. 

With these ideas of solar plasma flows and magnetic field dynamics in the background, 
I will trace the history of development of ideas in solar dynamo theory in the next section. 

2. Solar Dynamo Theory: Development of Ideas 

The behavior of magnetic fields in plasma systems (such as that in stellar interiors) is 
governed by the induction equation 

— = V X (v X B -r?V X B), (2.1) 
at 

where B is the magnetic field, v the velocity field and i] the effective magnetic diffusivity of 
the system - which in this case, encompasses the enhanced contribution due to turbulence 
in stelar interiors. Astrophysical plasma systems have a high characteristic magnetic 
Reynolds number (the ratio of the first to the second term on the R.H.S. of the above 
equation). The concept of flux-freezing holds in such systems (Alfven 1942), wherein, the 
magnetic fields remain frozen in the fluid flow - essentially because the diffusion timescale 
is much larger compared to the flow timescale. This allows the energy of convective flows 
in the solar convection zone (SCZ) to be drawn into producing and amplifying magnetic 
fields against dissipation - which is the dynamo mechanism mechanism in a nutshell! 

Assuming spherical symmetry (which is applicable to stellar interiors), the magnetic 
and velocity fields can be expressed as 

B = B^e^ + V X {Ae^) (2.2) 

V = r sin(6')f7e0 + Vp. (2-3) 

The first term on the R.H.S. of Eqn. 2.2 is the toroidal component (i.c, in the (/)-direction) 
and the second term is the poloidal component (i.e., in the r-9 plane) of the magnetic 
field. In the case of the velocity field (Eqn. 2.3), these two terms are the rotation and 
meridional circulation Vp, respectively. Due to the differentially rotating Sun, any pre- 
existing poloidal field would get stretched in the direction of rotation creating a toroidal 
component. It is believed that these strong toroidal flux tubes are stored and amplifled in 
the tachocline region - which coincides with a region of sub-adiabatic temperature gra- 
dient known as the overshoot layer and where magnetic buoyancy is suppressed (Spiegel 
& Weiss 1980; van Ballegooijen 1982). These horizontal toroidal flux tubes in the solar 
interior become unstable to magnetic buoyancy (Parker 1955a) when they come out in 
to the SCZ (through overshooting turbulence or meridional up-flows) and erupt through 
the solar surface creating bipolar sunspot pairs. During their buoyant rise, these flux 
tubes acquire tilt due to the Coriolis force generating the Joy's law distribution of tilt 
angles (D'Silva & Choudhuri 1993). 

For the dynamo chain of events to be to be completed, the toroidal held has to regen- 
erate the poloidal component so that the solar cycle can go on through a recycling of 
magnetic flux between the two components (in Eqn. 2.2) mediated via the plasma flows. 
The flrst mechanism proposed as a regeneration mechanism for the poloidal component 
of the magnetic field is the action of small-scale helical turbulence on rising magnetic 
flux tubes (Parker 1955b) - which twists the toroidal field back into the meridional place 
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Figure 3. A simulated solar butterfly diagram showing the evolution of the Sun's radial mag- 
netic field at the surface from a recently developed kinematic solar dynamo model. Eruptions 
of bipolar sunspot pairs are evident at lower latitudes; the redistribution of their fiux through 
surface flux transport processes involving meridional circulation, differential rotation and dif- 
fusion generates a large-scale global dipolar field. This cancels the older cycle polar fields at 
high latitudes and creates the polar fields of the new cycle. The observational features of the 
solar cycle are reproduced well through such simulations. This particular simulation is from a 
kinematic dynamo model that includes a realistic algorithm for the emergence of tilted bipolar 
sunspot pairs (Mufioz-Jaramillo, Nandy & Martens 2010). 

(thereby generating a r-9 or poloidal component). This mechanism relying on helical 
turbulence is traditionally known as the mean- field a-efFect. 

Simulations of the buoyant rise of toroidal flux tubes showed that the field strength 
of these flux tubes have to be close to 10^ G at the base of the SCZ to match the 
emergence properties of bipolar sunspot pairs (Choudhuri & Oilman 1987) and their tilt 
angle distribution (D'Silva & Choudhuri 1993; Fan, Fisher & DeLuca 1993). The required 
field strength is much stronger than the equipartition magnetic field strength in the SCZ 
(at which the energy in the fields and convection are in equipartition). This called into 
question the viability of the mean-field a-effect which is likely to get quenched by such 
strong fields. 

An alternative mechanism, originally proposed by Babcock (1961) and Leighton (1969), 
is now the leading candidate for the regeneration of the poloidal component of the mag- 
netic field. In this mechanism, the decay of tilted bipolar sunspot pairs, mediated via 
diffusion and meridional circulation creates a large-scale solar dipole from the emergence 
of individual bi-poles with a non-zero tilt (which have a net dipole moment). Thus the 
large-scale solar poloidal field is created and manifests as the global solar dipole config- 
uration - most intense at solar minimum - from which the next cycle toroidal field is 
generated. This keeps the cycle going. This observed process of surface flux distribution 
is reproduced by surface flux transport simulations (Mackay & Lockwood 2002; Jiang 
et al. 2011) and dynamo models (Munoz-Jaramillo, Nandy & Martens 2010) based on 
the Babcock and Leighton (BL) idea. More importantly, recently inferred correlations 
between observed proxies of the BL source term and the strength of the next sunspot 
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cycle (Dasi-Espuig et al. 2010) lend strong support to the BL mechanism for poloidal 
field creation. 

While full magnetohydrodynamic (MHD) simulations of the solar interior are becoming 
more and more sophisticated due to rapid improvements in our numerical capabilities 
they are still not capable of reproducing the rich variety of solar cycle observations. 
Kinematic dynamo models play a useful role here, wherein, their simplicity allows for 
a transparent underlying physics, and the capability to explore the full breadth of so- 
lar cycle properties through numerical investigations. In these models, the plasma flows 
are prescribed (and is not a dynamical variable) and the magnetic induction equation 
(Eqn. 2.1) is solved with suitable boundary conditions and parameterizations and al- 
gorithms necessary for accounting for the physics of the solar interior. Kinematic flux 
transport dynamo models (as opposed to models based on dynamo waves) rely on the 
plasma flows to transport flux between different source layers for the field components 
and have now become the leading tool for exploring solar activity (Jouve et al. 2008; see 
also the reviews by Charbonneau 2010) and recently developed models are good enough 
to reproduce not only the global features of the solar cycle, but also the dynamics of 
surface flux transport and polar field reversal accurately (Fig. 3). 

We do not provide a pedagogical description of the mean field dynamo equations 
here and instead refer (interested) readers to Nandy (2010a) - where a more complete 
description is provided, and to Nandy (2010b) - where some issues of current interest 
are discussed that will complement this review. I now jump ahead and devote the next 
section to discussing what I deem to be fertile grounds for dynamo theory in the coming 
years. 

3. A Future Outlook 

It is interesting times for solar dynamo theory. In the last decade major improve- 
ments have happened in the field of kinematic dynamo modeling and the so-called "flux 
transport" (kinematic) dynamo models have become useful tools for explaining various 
features of the solar cycle. / must note here in passing that I do not necessarily endorse 
the usage of the term "flux-transport" in this context as many dynamo modelers who de- 
clare their models to be "flux-transport" dynamo models do not include either magnetic 
buoyancy or downward pumping of magnetic flux - which are plausibly the most efficient 
of all flux transport mechanisms in the solar interior. Simultaneously, in the last several 
years or so, full MHD models have developed to the point where structured large-scale 
magnetic fields (Brown et al. 2011) and cyclic reversals in solar-like conditions are starting 
to be seen (Ghizaru, Charbonneau & Smolarkiewicz 2010). With the available long-term 
sunspot record in conjunction with polar magnetic field measurements stretching over 
three cycles and helioseismic observations covering at least a full cycle, we also now have 
the requisite observational data to constrain and indeed, drive dynamo simulations. I 
have no doubt that the most important advances in dynamo theory are going to come 
at the confluence of all of these, where kinematic and full MHD simulations of the solar 
cycle complement each other and assimilate magnetic field and plasma flow data. Below, 
I elaborate on some specific aspects - which are important and which I believe are going 
to garner sustained attention in the near future. 

3.1. It takes Two to Tango: Kinematic Dynamos and full MHD Simulations 
The history of solar dynamo theory has been such that kinematic dynamo modeling 
and full MHD simulations of solar and stellar interiors have often progressed in parallel 
with very limited feedback from each other. This is changing. An prominent example 
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is explorations of the nature of the meridional flow which is largely unconstrained by 
observations. Kinematic dynamo models, especially those based on the BL mechanism - 
wherein the source layers for the poloidal and toroidal field are spatially separated, find 
that the meridional circulation is an important factor in determining various aspects of 
the solar cycle, including fluctuations in the period (Charbonneau & Dikpati 2000) and 
the latitudinal distribution of sunspots (Nandy & Choudhuri 2002); the latter proposed 
that the flow has to penetrate below the base of the SCZ for reproducing the observed 
solar butterfly diagram. While this suggestion was controversial at the time, full MHD 
models are not beginning to find predominantly unicellular flows (averaged over solar cy- 
cle timescales) such as that used in kinematic dynamo models with overshooting plumes 
and plasma flows that can indeed penetrate below the base of the SCZ (Rogers, Glatz- 
maier & Jones 2006; Garaud & Brummell 2008). This penetration of plasma plumes and 
flows into the radiative interior is important because this region is stable to magnetic 
buoyancy; this in conjunction with the longer diffusion timescales there, allow for mag- 
netic field storage and slow transport (by the equatorward meridional counterfiow) that 
are relevant for establishing an 11 year solar cycle. Clearly, this convergence of kinematic 
dynamo models with MHD simulations of the solar interior reiterate the importance of 
meridional circulation vis-a-vis the solar cycle. 

It is to be noted that kinematic dynamo models does not necessarily self-consistently 
include the full range of physics possible in stellar interiors and therefore must assimilate 
those physical processes that full MHD simulations point out to be important. A case in 
point is the downward pumping of magnetic flux in the presence of rotating, stratified 
convection - a process that is very efficient at pumping magnetic fields down into the 
stable interior below the base of the SCZ (Tobias et al. 2001; Kapyla ct al. 2006). Only few 
kinematic dynamo models include downward pumping - those that do show that it has a 
significant influence on the dynamo (Guerrero & de Gouveia Dal Pino 2008); evidently, 
given its low timescales, the process of flux pumping competes with meridional flow and 
turbulent difi'usion in transporting flux between the source layers of the dynamo and 
should be considered in kinematic dynamo models. I believe useful exchanges between 
kinematic dynamo modeling and full MHD simiilations of the kind highlighted above are 
going to lead to more meaningful advances in the understanding of the solar cycle. 

3.2. Beyond Normal: Irregularities in the Solar Cycle 

Fig. 1 shows that the amplitude of the solar cycle varies significantly. While magnetic 
buoyancy presumably sets the order of magnitude field strength (Nandy 2002), smaller 
fiuctuations around these values likely have a different origin. Then there are grand min- 
ima phases like the Maunder minimum when the sunspot cycle basically stops. There 
are also variations in the period of the cycle, the length of the minimum between solar 
cycles and the strength of the polar field. Given that dynamo models are now reasonably 
successful in reproducing and explaining the regular features of the solar cycle, the chal- 
lenge ahead lies in confronting models of magnetic field generation with irregularities in 
the solar cycle. 

What causes these variations - is it due to stochastic fluctuations, e.g., in the meridional 
flow and the dynamo a-effect (Charbonneau & Dikpati 2000), or is it the non- linear 
feedback of the flelds on the flows (Tobias 1997) that may lead to chaotic modulation or 
is it due to sudden changes in some dynamo ingredient like the meridional flow (Karak 
2010). Some observational analysis seems to rule out the existence of low- dimensional 
chaos in the solar cycle (Mininni, Gomez & Mindlin 2002), however, the underlying non- 
linearities of the system can certainly exhibit chaotic modulation when the forcing is 
strong; the jury is still out on this one. 
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While it is tempting to surmise that the same process(es) that produce cycle to cycle 
fluctuations also sometimes bring down the activity levels to produce grand minima like 
episodes - that would be intellectually satisfying - does it hold in reality? The unusually 
long minimum of solar cycle 23 and a dynamo based explanation of this episode (Nandy, 
Munoz-Jaramillo, & Martens 2011) suggests that such variations in the nature of solar 
minima are caused by different processes than those that cause grand minima episodes. 
The way the Sun gets into Maunder-minima like episodes and the way it gets out from 
these quiescent phases also holds keys to understanding many of the subtleties of the 
dynamo process. We arc currently just scratching the surface here and need to know 
more about the origin and implications of solar cycle fluctuations and extreme events. 

3.3. Data Assimilation and Solar Cycle Forecasting 

The holy grail of solar dynamo theory is to be able to forecast the solar cycle. It is 
easier said than done, especially given the uncertainties regarding the relative roles of 

flux transport processes such as diffusion, circulation and downward flux pumping in the 
solar interior. Dynamo based predictions for the next cycle do not converge (Dikpati, de 
Toma & Oilman 2006; Choudhuri, Chatterjee & Jiang 2007) and this is due to differing 
importance laid on meridional circulation and turbulent diffusion (see Yeates, Nandy & 
Mackay 2008). The memory of the solar cycle is related to the timescale of the dominant 
flux transport process and since the build up of memory makes prediction possible, 
accurate predictions will require understanding the subtleties of various flux transport 
processes in the solar interior. Prediction also requires the assimilation of data. At the 
least, the data could be polar field (or surface magnetic field) observations - which 
are precursors of the next cycle (Dasi-Espuig et al. 2010) ~ and which can be fed into 
kinematic dynamo models to generate the toroidal field of the next cycle. 

However, more sophisticated treatment of the prediction problem, including predictions 
of solar cycle timing would also necessitate the realistic assimilation of observed plasma 
flow data in solar dynamo models; such efforts are just beginning (see Munoz-Jaramillo, 
Nandy & Martens 2009). Efforts are also on to take methods used in meteorology and 
develop them for solar cycle predictions (see e.g., Katiashvili & Kosovichev 2008 and 
Jouve, Brun & Talagrand 2011). Given the success of such techniques in weather pre- 
diction, this would seem to a pragmatic approach. In this approach, an observed state 
of the Sun is integrated forward in time using a solar model, and continuous corrections 
based on real-time observations are made to keep the model on track. Evidently, such an 
approach requires a good dynamo model to begin with and techniques for assimilating 
both plasma flow and magnetic fleld observations within this model. Givcm the need of 
reliable space weather forecasts, it would seem that attempts to predict the solar cycle 
would also drive the motivation to understand it better in the coming years! 

3.4. Solar-Stellar Connections: Stellar Activity and Dynamo Theory 
The Sun is one amongst many stars with similar properties. There is no reason than that 
a dynamo model for magnetic field generation would be unique to the Sun; in fact, one 
would think that the Sun just offers one data point for the dynamo, whilst observations 
of other solar- like stars offer multiple data points in the parameter space of e.g., differ- 
ential rotation or convective turn-over time - to constrain the dynamo mechanism with. 
An extreme form of this argument would imply that an acceptable dynamo model for 
the Sun should also be able to explain the activity of other stars with suitably adjusted 
parameters conforming to the interior of these stars. Typically however, not all dynamo 
models are confronted with stellar observational data. This is because the stellar data 
does not have well-constrained plasma flow observations (differential rotation and merid- 
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ional circulation) and in the absence of these, there are too many free parameters in 
spatially extended dynamo models for the whole exercise to be meaningful. Often, low- 
order mathematical models have been used instead to explore the wide range of stellar 
activity observations and they provide useful insight on how the exhibited dynamics is 
related to the underlying structure of the dynamo equations (see e.g., Wilmot-Smith et 
al. 2005, 2006 and references therein). 

Instead of modeling individual, unresolved stars, an alternative and promising ap- 
proach is to extend solar dynamo models to explore well-known relationships that exist 
in the stellar data, e.g., the rotation rate, activity amplitude and period relationship 
(Noyes, Vaughan & Weiss 1984; Brandenburg, Saar & Turpin 1998; Nandy 2004; Jouve, 
Brown & Brun 2010). Surface flux transport ideas gleaned from the Sun can also be 
adapted to explain the formation of polar caps and starspot dynamics (Schrijver & Title 
2001; Holzwarth, Mackay & Jardine 2006; Isik, Schiissler & Solanki 2007). Such ap- 
proaches sometimes throw up surprises; it appears that the BL dynamo as envisaged for 
the Sun runs into some difficulty when confronted with the activity-period relationship 
from stellar observations (Jouve, Brown & Brun 2010). Given that the BL dynamo is sup- 
ported both by theoretical models and solar observations are we missing something when 
we confront such models with stellar data - perhaps a realistic handling of downward 
flux pumping? A combination of full MHD simulations of stellar interiors and kinematic 
dynamo modeling would be required to address such vexing issues; certainly more effort 
should be devoted to using stellar activity relationships to constrain dynamo theory more 
meaningfully. 

3.5. Space Climate: The Solar Dynamo in Time 

I would like to end by bringing up a topic that has not received much attention but 
which is going to bcicome important in the years to come. This pertains to uncovering 
the long-term evolution of the solar dynamo over the life-time of the Sun and using this 
in conjunction with ideas of solar luminosity variations with radius and age to come up 
with a complete profile of the radiative history of the Sun. Theoretical and observational 
studies of solar-like stars at different main-sequence ages are going to play an important 
role in this regard (Nandy & Martens 2007). With increasing importance of understand- 
ing the causes and (;onsequences of global climate change, the need to figure out how 
our parent star has shaped our atmosphere and climate over time will also become an 
important quest. 

/ would like to acknowledge the Department of Science and Technology of the Gov- 
ernment of India for supporting my research through the Ramanujan Fellowship. I am 
grateful to the Argentinians - specifically for being wonderful hosts, and in general - for 
giving tango to the world. 
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